The Anatomy of Abelian and Non-Abelian Fractional Quantum Hall States 
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We deduce a new set of symmetries and relations between the coefBcients of the expansion of 
Abelian and Non-Abelian Fractional Quantum Hall (FQH) states in free (bosonic or fermionic) 
many-body states. Our rules allow to build an approximation of a FQEI model state with an overlap 
increasing with growing system size (that may sometimes reach unity!) while using a fraction of 
the original Hilbert space. We prove these symmetries by deriving a previously unknown recursion 
formula for all the coefficients of the Slater expansion of the Laughlin, Read Rezayi and many other 
states (all Jacks multiplied by Vandermonde determinants), which completely removes the current 
need for diagonalization procedures. 
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Model wavefunctions, such as the Laughlin Moore- 
Read (MR) 0, and Read-Rezayi (RR) 3 

states, have 

so far provided the key to understand the physics of 
Abelian and non-Abelian FQH phases. Despite their ex- 
plicit availability (in some cases) in terms of the electron 
positions, their expansions in the non- interacting basis of 
occupation number states (Slater determinants or mono- 
mials) is unknown and is considered intractable. As a 
result, quantitative studies of these states have heavily 
relied on exact diagonalization methods 

In usual numerical methods, one starts from a model 
[3| Hamiltonian, generates the Lowest Landau Level 
(LLL) Hilbert space and diagonalizes the Hamiltonian 
within this space. However, one immediately hits an in- 
surmountable barrier: diagonalizing a Hamiltonian ma- 
trix of a factorially growing size. The only known sym- 
metry of the coefficients of the expansion of a FQH state 
in Slater determinants is that of the angular momentum 
Lz —Lz'- on the sphere this reflects the indistinguisha- 
bility of the North and South poles. This symmetry, only 
valid when the (useless for quasihole excitations), 

roughly halves the size of the Hilbert space needed to con- 
struct the model FQH state. It would hence be very bene- 
ficial to discover other symmetry rules that the expansion 
coefficients satisfy. There exist several past attempts at 
identifying the coefficients of the free many body states in 
the simplest interacting Laughlin 1/3 state £,]. How- 
ever, those attempts can only obtain a small number 
0(I/iV!) of these coefficients, a statement which becomes 
painfully obvious when the maximum size Laughlin state 
generated by these methods is still much less than what 
can be achieved by exact diagonalization studies. 

In the present letter we discover a new series of rules 
that the coefficients of model FQH states satisfy. Our 
rules are valid for all bosonic and fermionic Jack poly- 
nomial states which include the Laughlin, MR, and RR 
series. We first obtain an explicit recurrence relation 
for the coefficients of a fermionic model FQH state that 
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FIG. 1: Example of the product rule (see text) for fermionic 
and bosonic MR States. The rule is similar in idea to the 
computation of Feynman disconnected diagrams. 



completely removes the need to use diagonalization tech- 
niques, hence removing the most important obstacle to 
achieving larger size model wavefunctions. We then use 
this to obtain several rules that the coefficients of the 
expansion satisfy. Some of these symmetries are reminis- 
cent of product rules in Feynman disconnected diagrams. 
One specific case of our rules is equivalent to the existence 
of perfect off-diagonal long range order (ODLRO); other 
cases explain small parts of the observed entanglement 
spectrum (ES) 0] ■ For Laughlin and MR, the application 
of these rules roughly halves the Hilbert space needed to 
produce these states. Still the overlap with the exact 
state is increasing with system size and can approach 
unity in the thermodynamic limit. We implement our 
new methods to present a proof of principle break of the 
current size barrier and obtain the MR state for 24 parti- 
cles, with Hilbert space « 20"* times larger than current 
exact diagonalization techniques. A future paper will use 
the new muscle to compute quasihole propagators and 
topological properties. 

All the FQH states discussed in this paper are squeezed 
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polynomials: in the development onto occupation num- 
ber (Fock basis) states (monomials for bosons and slaters 
for fermions), we only find orbital occupations obtained 
by a so-called squeezing operation on the root occupation 
configuration Q Q . We represent an angular momentum 
partition A with length £x < N as a, occupation-number 
configuration n(A) — {n„i(A), m = 0, 1, 2, . . .} of each of 



the LLL orbitals (/i™(z) = (27rm!2™)-i/2z™ exp(- 
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with angular momentum Lz = mh, where, for m > 0, 
f^m(A) is the multiplicity of m in A. It can be only or 1 
for fermionic FQH states, but can be any positive integer 
for bosonic ones. The formalism used here is thoroughly 
presented in [1, . It has been showed [1, 0| that the TV- 
particle bosonic Read-Rezayi series of states (which in- 
clude Laughlin and MR) are the r — 2 single Jack poly- 
nomials (Jacks) J^{zi, ...,Zn) of parameter a — — 
and root partition n(A) = [kO^~^k...kO^~^k]. Quasihole 
states can be written as coherent state superpositions 
of Jacks This identification proves their squeez- 

ing property and allows; unfortunately, since the Jacks 
are symmetric polynomials, only bosonic FQH states can 
be accessed through this procedure, while the fermionic 
states are the physically relevant ones. The problem 
is more severe: the Jacks are expanded in quantum- 
mechanically orthogonal free-boson many-body states 
(monomials) mx{zi, zn) = )/ Hm 

where Per is the permanent. When multiplied by a Van- 



dermonde determinant, the monomial does not become a 
Slater determinant, this being the property of the Schur 
functions. As such, to obtain the Slater decomposition of 
the fermionic FQH state J"(zi, zn) Y[i<:ji^i ~ ^j) oii^ 
would first have to transform the Jack from monomial to 
Schur basis. This transformation involves knowledge of 
all the Kostka numbers, a long-standing unsolved math- 
ematical problem with no known efficient algorithm. We 
depart from this hopeless route and attack the problem 
differently: we want to find the coefficients of the 
decomposition: 

N 

Sx{zi...zn) = J^B{zi...ZN)Y\_{zi-Zj) = ^ bxf,s\^ (1) 

i<j fJ-^^ 

sU = Detiz:;'') is the quantum mechanically (unnormal- 
ized) orthogonal Slater determinant; all partitions fi are 
squeezed from the squeezed partition A 1, related to 
the Jack bosonic partition A^ by: Ai = Af -I- (TV — i). 
By rescaling we choose 5a a ~ 1- J^b is an eigcn- 
state of the Laplace-Beltrami operator HLB{a) with 
known eigenvalue [llj. It is then only a matter of (te- 
dious) algebra to find the new operator that diagonal- 
izes S'"(zi...Z7v). We call it Fermionic Laplace-Beltrami 
H[g{a)Sx{zi...ZN) = Ex{a)Sxizi...ZN), with Ex = 
E,A,(A,-2(i-l)z) + (i-l)((7V + l)E,A.-TV(iV-l)): 



We now consider the action of this operator on a Slater 
sl^i . The action of the kinetic term is J2i ^kslp — 
J2i Aifsl^. The action of the remaining "interaction" part 
is complicated. We sketch the significant step needed to 
untangle it. As the "interaction" is a two-body operator, 
it is revealing to consider its action on the two particle 
Slater sl^=(^^^^2) = ^i^z!^^ — Z2^Zi^- Its action on larger 
size Slaters can be decomposed in pairwise actions (with 
sign-keeping complications). We find i/isl(^j_^2) equals: 

2 

(Mi -M2 -2)sl(^i^^,) -)-2 ^ [ni- ^2-2s)s\i,^_s,t,2+s) 

s=l 

The interaction term acting on a Slater creates all pair- 
wise squeezed configurations of a partition, but its ac- 
tion is different than that of the corresponding Laplace- 
Beltrami operator on a monomials. When generalized 
to any number of particles, by virtue of dealing with 
fermions, the above action acquires a —1 sign every time 
that either fii — s or fij + s crosses another partition com- 




ponent /ifc. Taking the action of H^g{a) on S*" in Eq(IT]), 
we rescale the summation over s above and use linear 
independence of Slaters to obtain: 

The sum is over all partition 9 = {fii, ...,fii + s, — 
s,...,/iAr) which strictly dominate = (/ii, ^tv) and 
are squeezed from the root A. Ngw is equal to the num- 
ber of swaps required to order the partition 6. Px{ct) = 
X^i Ai(Ai + 2i{l — 1/a)). Eq.([3]) differs from a known re- 
lation for the exp ansion of a Jack polynomial in terms 
of monomials besides different Pa (a), the factor 

{^i — fij) in the summand of Eq([3]) does not depend on 
the partition 6 (does not depend on s) whereas it does in 
the bosonic case {^i — fij + 2s). For a = —{k + 1), Eq.Q 
gives the coefficients of the RR fermionic states. It com- 
pletely removes the need for diagonalization, coefficients 
being computed iteratively starting from the root par- 
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tition The algorithm's numerical stability is high 

and has been checked through some coefHcients which 
are explicitly known. Largest system sizes previously ob- 
tained using diagonalization can now be done on a single 
CPU workstation in a matter of a few hours (for the 1/3 
Laughlin state at = 15, the new algorithm requires 
500 times less CPU time). The main bottleneck is now 
the Hilbert space storage which may still be huge. 

We now give an example and sketch the proof of a sim- 
ple rule we found the coefficients of a model FQH state 
to satisfy. As above, we work in the un-normalized or- 
bital basis of Slaters or monomials; this completely sepa- 
rates the polynomial part of the problem from the scalar- 
product part, which involves geometry-dependent nor- 
malizations and is hence not fundamental. Once the state 
is obtained in un-normalized basis, it can be trivially 
transformed in normalized basis depending on whether 
we want to perform calculations on the sphere, disk, 
cylinder or any other genus geometry. 

The Moore-Read state for A*" = 10 particles is a lin- 
ear combination of Slater determinants squeezed from 
the root configuration n(A) = 110011001100110011 
[1, 0]. Let us determine the coefficient of e.g 
101101000101111010 (see FigP). This configuration has 
the special property that, if cut in two after the first 7 
orbitals, each resulting partition can be squeezed from 
its own root partition: at the left of the cut, 1011010 
can obtained by squeezing on the root partition = 4 
MR root occupation 1100110; at the right of the cut, 
00101111010 can be squeezed from the = 6 MR root 
occupation 01100110011 (extra fluxes -zeroes- can be ne- 
glected). We proved a product rule: the coefficient of 
101101000101111010 in the A^ = 10 particle MR state 
(140) equals the product of the coefficients of the two 
disconnected A^ = 4 and A^ = 6 pieces of the partition: 
—2 and —70. A similar property occurs in the bosonic 
MR state (see Fig[Tj), and for any disconnected occupa- 
tion numbers. The rule also allows further "dissection" of 
FQH states. It can be trivially generalized to any product 
of disconnected pieces. This is a direct consequence of the 
recursive property of the product rule : some coefficients 
might be computed from those of two smaller systems 
which in some case, could be deduced from those of two 
even smaller systems, etc. This rule bears resemblance 
to the rules of computation of Feynman disconnected di- 
agrams, where the value of a diagram that can be cut in 
two disconnected parts can be obtained as the product 
of the values of the two disconnected parts. 

This simple rule valid for all Jack polynomials, bosonic 
and fermionic, at any a, of any partition. We now very 
roughly sketch the proof for the fermionic Jacks using our 
newly found formula. For bosonic Jacks, a similar proof 
can be obtained. Assume we want to determine the coeffi- 
cient, in an arbitrary fermionic polynomial S'^{zi, zn), 
of a configuration ^ < A that can be divided in two dis- 
connected sets. By assumption fi = {ijLA,y^B) of Na and 



Nb particles [Na + Nb = A"), is supposed to be divis- 
ible in two independent partitions: ^a, squeezed from 
the root partition of Na particles A^i = (Ai, A^r^), 
and /is, squeezed from the root partition of Nb particles 
As = (Aat^+i, Aat). The proof proceeds by induction. 
Assume the product rule is valid for ALL partitions, of 
"separable" form 9 = {9a, 9b) where ^ia < 9a < and 
fJ.B < 9b < ^B- The coefficients of 5*^ are given in EqU]) 
as a recursion of coefficients from partitions which dom- 
inate fi. Crucially, if /i has the separable form chosen 
by us, then the partitions which dominate fi and which 
enter EqlT]) are also separable, of the form 9 = {9a, 9b)- 
Hence the sum in Eq((T]) separates in two distinct sums of 
(/ii— /Xj)&;je(a)-(— l)'^™ over the disconnected parts of the 

partitions 9: Ee-a.<e<x = E,,^<e^<A^ + E,,^<e^<A« ■ 
In the first sum /i^ , fj,j belong only to the left hand side 
9a of the partition and not to hb- the partition in the 
first sum reads {9a, A^b)- In the second sum they belong 
only to the right hand side 9b of the partition, not to 
/i^:the partition in the second sum reads {^a,9b)- By 
the induction assumption, all partitions which dominate 
satisfy the product rule: bx{0^^^s){a) = foA(e^,AB)(a) ■ 
bx(Xji^fiB){(x) where {9a, As) is denoted as being the par- 
tition formed by 9a and the ground-state partition As, 
and similarly for the other coefficient. The sums then 
become E^^<0^<a^ (A^i " Mj>A(e^,MB)(a) ' (-1)^™ = 

&A(AA,MB)(a) • E^^<fl^<A^(M« - Aij)&A(e^,AB)(a) ' (-1)^™ 
and similar for the sum over 9b- We notice that the sum 
is proportional to bx[^^,\g)- By using: 

PA(a) = PA(a) -pLab(") -Paamb(") 

we can prove that bxt,{a) = bx^^^^XsM) ' ^a(Aa,pb) 
which is the product rule described. It is valid for all 
Jacks, at any a and hence for all Read Rezayi states. We 
conjecture it is also valid for some non-Jack states [Hj. 

In special cases, wc can identify the product rule with 
physical properties of the FQH states. Cut the A^-particle 
Laughlin v = \/2 state 101010101... 101 into two pieces, 
one of 2 orbitals and 1 particle (10) and the second of 
N - I particles 1010101. ..101. The product rule shows 
that the coefficients of any configuration in the A^-particle 
state that has 1 particle in the first orbital has the same 
coefficient as the one that the configuration obtained by 
deleting 10 from the A^-particle state would have in the 
A^— 1 particle state. Since deleting 10 from the Af-particle 
state is equivalent to applying the operator ft,^(z)i/)(z) 
where h{z) — {z — ZiY , in this specific case the 
product rule can be re-written ji/'Ar-i) = {z)iP{z)\iI)n) , 
which is ODLRO statement. In the non-abelian case, a 
similar statement occurs, when the ODLRO concept is 
generalized to non-abelian states [lS|. 

The product rule also explains aspects of the ES of the 
FQH states. On the sphere, we cut the state into two 
hemisphere blocks A and B- Following 0], we introduce 
the ES ^ as Ai = ea;p(— ^i), where \i are the eigenvalues 
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of the reduced density matrix pA of one hemisphere. The 
eigenvalues can be classified by the number of fermions 
Na in the A block, and also by the total angular momen- 
tum lI"^-* of the A block. It was empirically found for the 



MR01, Laughhn [16|, and Jain Composite Fermion [17 1 
states that the low-lying spectrum of the reduced den- 
sity matrix for fixed Na, plotted as a function of , 
displays a structure refiecting the edge CFT. In particu- 
lar, at maximum L^'^'' one finds one single eigenvalue of 
the ES, irrespective of where in orbital space the cut was 
performed. A generic state is expected to have an order 
min(Af^!, [N — Na)^-) levels. Our product rule proves that 
only one eigenvalue is possible at the maximum Li^). As 

(A) 

we consider the maximum Lz , the density matrix at 
this angular momentum accesses states that are squeezed 
separately within each hemisphere; if a state would have 
a configurations squeezed across the cut, then the angular 
momentum ii"^^ would no longer be maximum. As the 
states are squeezed separately within each hemisphere, 
the product rule renders the rows and columns of the 
density matrix linearly dependent, and hence 1 non-zero 
eigenvalue of the Schmidt spectrum. 

Given its generality, one can ask whether the product 
rule works for any L = state in the LLL. It does not. 
A counter example is the linear superposition of the two 
states obtained by applying the L = Q condition on the 
subspace of states squeezed from the Jain state root par- 
tition 2010110102 17]. Cutting the state symmetrically 
in half after 5 orbitals, the ES generically has 2 eigen- 
values at the largest lI"^' = 12, and the product rule 
in general doesn't apply. However, since the counting in 
the ES is conjectured to be equivalent to the counting of 
edge-modes, and since the Jain state has 1 excitation at 
maximum L^^^ - which we showed is equivalent to the 
existence of a product rule - the rule should work for the 
Jain state as well in the thermodynamic limit. 

How much of a state can be obtained by the prod- 
uct rule? Assume that the N particle Laughlin state 
"^10101 ioioi(-^i' •^^^) known (this trivially implies 
that the Laughlin states for smaller number of par- 
ticles are known). In the + 1 Laughlin state 



1010101. ..1010 

not be obtained using the product rule are dominated 
by the configuration where the 2 particles in the North 
Pole and South Pole orbital are squeezed together once: 
[01 10101... 10101 10]. Hence, the number of configura- 
tions that one can obtain in the 1/2 Laughlin state is 
roughly 45-50% in the thermodynamic limit. For non- 
abelian states, we can obtain a smaller percentage - for 
MR, we get roughly 30% in the thermodynamic limit. 
We then compute the overlap of the exact iV-particlc 
Laughlin state with the state obtained by applying our 
product rule on the exact iV — 1 Laughlin state. We ob- 
tain 0.9875 for = 6 up to 0.9977 for N = 15; note that 
the overlap grows, and likely reaches 1 in the thermody- 



namic limit, even though we miss half the configurations 
of the FQH state. Similarly for the Moore- Read state, the 
product rule cannot obtain configurations squeezed from 
[112020202. ..2020211]. The overlaps are 0.8858 for A^ = 8 
up to 0.9383 for A^ — 22. For fermionic states, the over- 
laps are a bit lower (like 0.9603 for the 1/3 Laughlin state 
at A^ = 15) but with smaller fraction of the Hilbert space 
and still increase with system size. We have checked that 
the product state entanglement entropy is less than 5% 
higher compared to the model state. 

In fact, we can do better. We have also proved a series 
of rules for configurations that cannot be obtained by 
the product rule. For example, for Laughlin v = 1/2, the 
product rule is useless in obtaining the coefficients of any 
of the configurations squeezed from 01101010. ..1010110. 
However, we have proved another rule that can obtain 
some of these coefficients. We can prove that the coeffi- 
cient of the partition OllOA in the A^ + 1 particle state 
(where A is any partition squeezed from 1010. ..1010110) 
is proportional to the coefficient of the partition 01 A in 
the A^ particle state. The proportionality constant can 
be obtained by considering the first two members of these 
partitions 0110101. ...010110 and 100110101. ..010110. For 
Laughlin 1/2 state, the proportionality constant is — 1. 
The same rule works for non-abelian states. For example, 
for the MR state, the product rule is useless in obtain- 
ing configurations squeezed from 1120202. ..2020211; our 
new rule proves that the coefficient of the configuration 
11 20 A (where A is squeezed from 202...2020211) in the 
A^ -f 2-particle state is proportional to the coefficient of 
the configuration llA in the A^-particle state. For MR 
A — 22, using this rule in addition to the product rules 
improves the overlap to 0.9417 while adding only 1.2% 
of the vector components. The rule above is part of a 
unified rule that allows the determination of coefficients 
of all configurations where the particles in the North and 
South pole have been squeezed only once; using the two 
rules, for Laughlin v — 1/2 we can obtain all the coeffi- 
cients for all configurations except the ones squeezed from 
0020101. .1010200; for MR, we can obtain all coefficients 
except the ones squeezed from 0220202... 02020220. The 
new rule explains the existence of only 1 eigenvalue in 



j^(zi, ...,zjv+i) the configurations that can- the ES at L 



(A) _ 



1, and, unlike the product rule. 



is only valid for FQH states and their quasiholes (not for 
all Jacks of arbitrary a). Due to lack of space, we will 
present this rule, as well as other conjectured rules in a 
separate long publication. 

Conclusion We have obtained a new formula for the co- 
efficients of the expansion of fermionic model FQH states 
in Slater determinants that completely removes the need 
for matrix diagonalization. Using this formula, we have 
proved a series of new rules of these coefficients that al- 
low further massive reductions in the Hilbert space size 
necessary to build a state. We hope the rules presented 
here will redefine the way numerics on model FQH states 
is currently done. 
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